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^ ■ ABSTRACT 

Q . We have conducted the first systematic study of Herbig Ae/Be stars using the tech- 

nique of long baseline stellar interferometry in the near-infrared, with the objective of 
j§ \ characterizing the distribution and properties of the circumstellar dust responsible for 

the excess near-infrared fluxes from these systems. The observations for this work have 
^ I been conducted at the Infrared Optical Telescope Array (IOTA). The principal result 

' of this paper is that the interferometer resolves the source of infrared excess in 11 of the 

^ 15 systems surveyed. A new binary, MWC 361-A, has been detected interferometrically 

for the first time. 

The visibility data for all the sources has been interpreted within the context of four 
simple models which represent a range of plausible representations for the brightness 
distribution of the source of excess emission: a Gaussian, a narrow uniform ring, a flat 
blackbody disk with a single temperature power law, and an infrared companion. We 
find that the characteristic sizes of the near-infrared emitting regions are larger than 
previously thought (0.5 — 5.9 AU, as given by the FWHM of the Gaussian intensity). 
A further major result of this paper is that the sizes measured, when combined with 
the observed spectral energy distributions, essentially rule out accretion disk models 
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represented by blackbody disks with the canonical T(r) oc r~^/^ law. We also find that, 
within the range observed in this study, none of the sources (except the new binary) 
shows varying visibilities as the orientation of the interferometer baseline changes. This 
is the expected behaviour for sources which appear circularly symmetric on the sky, and 
for the sources with the largest baseline position angle coverage (AB Aur, MWC 1080- A) 
asymmetric brightness distributions (such as inclined disks or binaries) become highly 
unlikely. 

Taken as an ensemble, with no clear evidence in favor of axi-symmetric structure, 
the observations favor the interpretation that the circumstellar dust is distributed in 
spherical envelopes (the Gaussian model) or thin shells (the ring model). This inter- 
pretation is also supported by the result that the measured sizes, combined with the 
excess near-infrared fluxes, imply emission of finite optical depth, as required by the 
fact that the central stars are optically visible. The measured sizes and brightnesses 
do not correlate strongly with the luminosity of the central star. Moreover, in two 
cases, the same excess is observed from circumstellar structures that differ in size by 
more than a factor of two, and surround essentially identical stars. Therefore, different 
physical mechanisms for the near-infrared emission may be at work in different cases, or 
alternatively, a single underlying mechanism with the property that the same infrared 
excess is produced on very different physical scales. 

Subject headings: instrumentation: detectors, interferometers — binaries: close — 
techniques: interferometric — circumstellar matter — stars: early type, emission line, 
formation, imaging, pre-main sequence — infrared: stars 

1. INTRODUCTION 

The Herbig Ae/Be (HAEBE) stars are understood to be young stellar objects of intermediate 
mass (1.5 < M/Mq < 10) (Herbig 1960). Observationally, they are defined by spectral types B to 
F8 with emission lines and the presence of infrared (IR) to sub-mm excess flux due to hot and cool 
dust. The pre-main sequence nature of HAEBE stars is now well established, based on comparison 
of their location in the HR diagram with theoretical evolutionary tracks (Strom et al. 1972; Cohen 
& Kuhi 1979; van den Ancker et al. 1998; Palla & Stahler 1991), their physical association with 
the dark clouds in which most HAEBE stars are found (Finkenzeller & Jankovics 1984), and their 
high rotational velocities (Finkenzeller 1985). The question of the geometrical distribution of the 
cool material which gives rise to the characteristic long wavelength excess, however, has remained 
controversial and is still the subject of considerable debate. 

In the case of the T Tauri stars, the excess continuum flux is successfully reproduced without 
extinguishing the central star at visible wavelengths by a model in which a circumstellar (CS) disk 
is heated both by the central star and by active accretion. The disk is also the source of the 



-3- 



observed activity (permitted and forbidden emission line intensities and assymetric profiles, UV 
excess, veiling of absorption lines) through magnetospheric infall onto the stellar surface (Muzerolle, 
Hartmann, & Calvet 1998). An important consideration in understanding the success of this model 
in this case is the fact that the level of activity observed can not be explained by any known 
mechanism involving the star alone. 

For the HAEBE stars on the other hand, the situation is considerably more confused, since 
the much higher stellar luminosities permit a number of alternative physical processes to reproduce 
the observations. Moreover, the much smaller number of sources and their more rapid evolution 
make statistical analyses less reliable. As an example, Bohm &; Catala (1995) conclude that the 
origin of outflowing wind in HAEBE stars arises in an accretion disk boundary layer as in the T 
Tauri stars, whereas Corcoran & Ray (1998) favor the interpretation that the outflow arises in the 
star itself. Likewise, the interpretation of the line profiles of forbidden lines is controversial. In the 
case of the classical TTS, nearly all forbidden lines lack a redshifted component, a fact attributed 
to occultation of the receding part of the wind by the opaque CS disk. A similar analysis for the 
HAEBE stars however shows that the forbidden lines are symmetric and therefore inconsistent with 
such a model (Bohm & Catala 1994), although Corcoran & Ray (1997) reach opposite conclusions 
based on similar data. 

In a much cited paper, Hillenbrand et al. (1992) (hereafter HSVK) classified 47 HAEBE stars in 
three groups according to the shape of their spectral energy distributions (SED). Group I contained 
stars with IR spectral distributions XFx ~ A~^/^, reminiscent of the spectral shape expected for a 
star - disk system in which the disk both reprocesses starlight and/or is self luminous via active 
accretion onto the central star (Basri &; Bertout 1993). Group II contained stars with flat or 
rising IR spectra, which require the presence of gas and dust not confined to a flat disk. Finally, 
Group III contained stars with small or no IR excess, similar to the classical Be stars, in which 
the small excesses above photospheric levels arises from free-free emission in a gaseous disk or 
envelope. HSVK modelled the SEDs of Group I stars assuming that the IR excess does in fact 
arise in a reprocessing and actively accreting disk, and they proposed an evolutionary sequence 
Group II Group I Group III, where the originally complex gas and dust environment that 
results from the star formation process leads natTirally to the formation of an accretion disk which 
gradually disappears. HSVK found that the SEDs of Group I sources could be satisfactorily fit 
provided that the accretion rates are relatively high (6 x 10^'' < M < 8 x 10^^ Mq yr^^) and that 
opacity holes of a few stellar radii exists in the inner region of the disk. These opacity holes were 
tentatively interpreted as being due to dust sublimation or, alternatively, to the interaction of the 
stellar magnetosphere with the disk, which lifts gas from the disk and forces accreted gas and dust 
to proceed along the stellar magnetic field lines (Konigl 1991). Although the same conclusions 
were reached by Hartmann et al. (1993), these authors also pointed out the fact that for the high 
accretion rates required to fit the SEDs, the gas in the region interior to the inner hole would 
be optically thick, resulting in additional emission which is not observed. Moreover, the derived 
accretion rates are inconsistent with the observed absence of optical veiling in the spectra of HAEBE 
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stars (Ghandour et al. 1994). The magnetosphere hypothesis, on the other hand, requires stellar 
magnetic fields strengths (~ 10 KG) not normally observed in early-type stars. 

At the very least then, it can be said that, although it would appear natural to assume that 
similar physical processes occur in the CS environment of T Tauri and HAEBE stars, the existence 
of accretion disks around HAEBE stars is far from being firmly established. As a consequence, 
alternative ideas have been considered to explain the observed IR excess, such as symmetric dust 
(and gas) envelopes (Berrilli et al. 1992; Pezzuto et al. 1997). It must be said that in this case too 
there is lack of consensus about the precise properties of these envelopes (density profiles, emission 
mechanisms), or even the question of whether or not envelopes are capable at all of explaining the 
observed SEDs (see for example Miroshnichenko et al. 1997 and Natta et al. 1993). One recurrent 
difficulty with symmetrical envelopes is that they tend to result in higher visual extinctions toward 
the central star than are observed, which has led some authors to postulate special geometries 
which provide a low optical depth line of sight to the star (Hartmann et al. 1993), or to models 
of more tenous envelopes consisting of very small, transiently heated dust grains Hartmann et al. 
(1993). 

In summary, there is currently no physically accepted model to account for the IR emission 
of HAEBE stars. At the same time, the answers to many questions on the evolutionary status 
and the origin of the activity and variability depend critically on the relative importance of CS 
distribution of material in disks or envelopes at different spatial scales. A desirable approach is 
to constrain the models not only with respect to the observed SEDs, but also with respect to the 
spatial distribution of the IR emission. Clearly, one way to make progress toward this important 
question is to obtain observations of the near-IR emission with sufficient angular resolution to 
resolve the structure around the central star, and that is the principal goal of the work reported 
here. 

In § 2 we describe our experimental procedure and data analysis method, including a novel 
technique for calibration of visibility amplitude data in the presence of atmospheric tip-tilt wave- 
front errors. In § 3 we present the list of HAEBE stars observed, log of observations and the 
calibrated visibility data. In § 4 we decompose the measured total fluxes into stellar and excess 
near-IR fluxes, an essential step before the visibility data for each source can be modelled in § 5. In 
§ 6 we summarize our results, and consider the ensemble of data in order to discuss the implications 
of our results for the nature of the CS environment in HAEBE stars. 



2. EXPERIMENTAL PROCEDURE AND DATA REDUCTION 

2.1. Experimental Procedure 

The observations described in this paper were carried out at the IOTA, a Michelson stellar 
interferometer located on Mount Hopkins, Arizona (see Traub 1998 for a recent review of the 
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lOTA instrument). Observations were made in the near-IR H (Aq = 1.65//m, AA = 0.30/xm) and 
K' (Aq = 2.16/um,AA = 0.32;um) bands; and using two IOTA baselines, of lengths B = 21 m 
(North-South orientation) and 38 m (North-North East orientation). For reference, the resolution 
corresponding to the longest baseline, as measured by the full-width at half maximum (FWHM) 
of the response to a point source, is \/2B = 4 mas (1.8 AU) and 6 mas (2.8 AU) at H and K' 
respectively, where the linear separations are given for a median distance of 460 pc to the stars in 
our sample. 

At the IOTA, starlight is collimated into 4.5 cm diameter beams by the telescope assemblies, 
of aperture diameter D = 45 cm, and are directed in vacuum towards the delay line carriages. The 
delay lines serve to maintain the zero optical path difference (OPD) condition between the two 
arms of the interferometer, from the source to the beam combination point, by reflecting the beam 
that requires a delay onto movable dihedral mirrors which track the sidereal motion of the target 
star. The path-equalized beams then emerge from the vaccuum tank and propagate in air inside 
the beam combination laboratory. In this experiment, a pair of dichroics reflect the near-IR light 
towards an optical table containing the beam combination optics and fringe detection system. The 
visible light in the beams is transmitted toward the CCD-based wavefront tip-tilt compensation 
servo system. 

Interference is produced in the pupil plane, by combining the collimated and path-equalized 
beams at a beam splitter. The two outputs of the beam splitter are focused onto two separate 
pixels of a NICM0S3 array, and a scan containing an intensity fringe packet is recorded at each 
pixel by modulating the OPD by ±60 fim via an extra reflection in one of the arms at a mirror 
mounted on a piezo stack driven by a highly linear triangle waveform. The scan rate is selectable 
in the range 1 — 10 scans/sec. 

Each scan contains 256 samples, and the integration time per data point is in the range 
0.18 — 2.0 ms, chosen as the optimum in order to achieve maximum flux sensitivity while keeping 
the scan time short in order to freeze the wavefront piston fluctuations induced by the atmosphere 
during the fringe packet acquisition. With this detection system, the limiting magnitudes for real- 
time visual detection of a fringe packet in a single scan are 7 at H-band (read-noise limited) and 
6.2 at K'-band (background limited), given here as the magnitudes of the faintest sources observed 
under typical conditions (Millan-Gabet et al. 1999a; Millan-Gabet 1999). 

A typical observation consists of 500 scans obtained in 1-8 min, followed by 10 scans on the 
sky which are used to measure the background flux. Target observations are interleaved with an 
identical sequence obtained on an unresolved star, which serves to calibrate the interferometer's 
instrumental response and the effect of atmospheric seeing on the visibility amplitudes. The target 
and calibrator sources are typically separated on the sky by a few degrees, and observed a few 
minutes apart; these conditions insure that the calibrator observations provide a good estimate of 
the instrument's transfer function. 
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2.2. 



DATA REDUCTION AND CALIBRATION 



2.2.1. Fringe Fitting 



The data reduction process begins by estimating the fringe visibihty, i.e. the fringe amphtude 
relative to the background-subtracted mean intensity in each interferogram. The raw response from 
the NICM0S3 pixels is first corrected for their intrinsic non-linearity, for both star and sky scans. 
A "median sky scan" is then computed as a point-by-point median of the 10 individual sky scans, 
and subtracted point-by-point from the star scans. If we call the resulting scans I a and Ib for each 
side of the beam splitter, we then compute reduced scans according to 



where / denotes the mean. The above operation: (1) compensates for any existing photometric 
unbalance between pixels A and B; (2) improves the SNR, since the beam splitter outputs Ia and 

Ib are out of phase by vr, while the noise is uncorrclatcd; and (3) eliminates common-mode noise, 
most notably intensity fiuctuations in the signal due to scintillation and residual seeing- induced 
motion of the focused star images outside the NICM0S3 pixels. 

The central 3 fringes in the reduced scans are then fitted in the time domain to a point 
source response template computed as the Fourier transform of the spectral bandpass used in the 
observation (H or K'). As a result of the fringe fitting procedure, we obtain the following quantities 
for each reduced scan: mean, fringe amplitude, position of the central fringe in the scan, and 
sampling rate. The data are flagged as being a false fringe identification if an individual fringe 
visibility is outside the range — 100%, or if either the central fringe position or sampling rate 
has a value which is more than three times the standard deviation obtained from the ensemble of 
the 500 reduced scans in the observation. Finally, an entire observation is rejected if the standard 
deviation of the central fringe position values exceeds L/2, where L = 120 /xm is the optical path 
corresponding to the scan length. 



Propagation through the Earth's atmosphere and inside the instrument degrades the inter- 
fering wavefronts differentially from their initial perfectly plane shapes, and therefore reduces the 
fringe visibility that is measured. The instrumental terms are constant or slowly changing and can 
therefore be monitored and calibrated using observations of a reference star known to be unresolved 
or of well known angular diameter. The atmospheric terms, however, change on time scales much 
shorter than the time required to do an observation and therefore cannot be calibrated using ob- 
servations of a reference star. In this section we address the problem of estimating the visibility for 
an observation in a way that is unbiased by low-order atmospheric effects. 



/ = 



Ia/Ia - Ib/Ib 
I a/ I A + Ib/Ib 



(1) 



2.2.2. Visibility Estimation 
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A histogram of visibility values obtained in a typical observation is shown in Figure 1, along 

with a curve representing the Gaussian distribution with the measured mean and variance. It can be 
clearly seen that the data distribution is not well represented by the Gaussian distribution, so that 
traditional estimators such as the mean would not be appropriate. For this work, we have developed 
a new technique which relies on our capability to record hundreds of individual fringes in only a 
few minutes, so that the distribution of visibilities is well determined and used in the estimation 
process. The size of the IOTA telescopes is well matched to the area over which an atmospherically 
disturbed wavcfront is expected to remain approximately flat, nominally ro ^ 50 cm in the near- 
IR. Therefore, the atmospheric visibility degradation is likely to be dominated by the residual 
differential wavefront tilts uncompensated for by the tip-tilt servo. Consequently, we consider the 
following representation for the measured visibilities 

^measured ~ ^true ' Vinst ' [^iU + V^oise] i = ^, ■ ■ ■ ,n (2) 

where Vtme is the true object visibility, Vinst represents the effect of slowly changing instrumental 
terms, such as wavcfront curvature or polarization mismatch, V^-^ is the visibility that results from 
interfering two tilted wavefronts, V^oj^e results from the fringe fitting process in the presence of 
read and photon noise, and n = 500 is the number of individual visibility values. Considering 
linear phase gradients on each aperture of differential slopes along orthogonal directions, it 

is straightforward to show that the resulting tip-tilt term is given by 

^ut= '^'^^l^f^ ' ^itl^^ A<^; = ^{AcPl + Acl^lY (3) 

A realization of the visibility distribution is generated by drawing values, n times, of 
and ^jioise from Gaussian distribiitions of [mcan,variance] equal to [Acpr , c^^^] and [0 , o'y^^_^J 
respectively. The physical motivation is that the mean differential tilt (A(/),.) represents a fixed 
instrumental misalignment, affected by seeing- induced fluctuations of variance <7%^^^- 

For a each observation, an estimate of the product V = Vtme ■ Vmst is obtained by finding the 
best match, in a sense, between the data distribution and the above model, in a grid generated 
for the complete range of the 4 parameters: 14r«e) '^y„oise' and a^^^. In practice, due to 
the fact that in this model it is difficult to distinguish the effects of the true visibility and mean 
tilt parameters on the distribution, the value of the mean tilt for a given night and for a given 
experimental set-up is obtained from the mean of preliminary fits to all the calibrator observations, 
and fixed to this value in subsequent fits of the target and calibrator observations. Similarly, the 
noise parameter is constrained for each observation to be the RMS in the visibilities obtained in 
the fringe fitting process. An example of the best-fit model distribution found is also shown in 
Figure 1. 
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Fig. 1. — Comparison of visibility data in an observation (500 scans) with a Gaussian distribution 
and the visibility distribution model (VDM) . In this case the visibility estimate given by the mean is 
42.8±0.6%, while the estimate given by the VDM is F = 60.9±0.3%. The other VDM parameters 
obtained are: Acf)^ = 0.15 rad, (Ja</,^ = 0.15 rad and (Tv„oise — 



2.2.3. Calibration 

A spline interpolation between sequential V estimates obtained in observations of calibrator 
sources, for which Vtme is known, yields a curve which reflects the time evolution of the instrument. 

Interleaved target observations, taken under identical instrumental conditions, are calibrated for 
the Vinst term by dividing the target estimate by the interpolated calibrator estimate: V^^^^'^^ = 
^true ' V^"'^^^^ /V'^"'^ The spectral types and visual magnitudes of all the calibrator stars used in 
this paper (with the possible exception of HR8881, see section 3) have angular diameters less than 
1 — 2 mas, therefore, we use V^^ue = 



2.2.4. Errors 

The internal uncertainty from fitting individual fringes is generally small (<1%), as is the 
formal uncertainty in the visibility for an observation obtained from the fit to the visibility distri- 
bution model (^ 2 — 3%). However, the RMS of independent visibility measurements on the same 
object and at similar spatial frequencies is larger, typically <5% p-p (pcak-to-pcak) within a single 
night of observation, although larger, <15% p-p in H-band data and <20% p-p in K'-band data, 
systematic differences occasionally exist when comparing data from different nights, both of which 
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indicate the presence of low-frequency calibration errors. A likely explanation for this excess error 
is that our expectation that rQ{H,K') ~ D does not strictly hold, and that wc see the effects of 
un-modelled wavefront curvature errors induced by the turbulent atmosphere. Therefore, the total 
uncertainty on an individual visibility measurement due to all sources is best represented by the 
RMS of independent measurements separated by hours of days. Using this measure, we find that 
the use of the visibility distribution model yields more consistent results, reducing the scatter of 
independent measurements by factors of 2 — 3 compared to the mean estimator. This improvement 
is illustrated in Figure 2, which compares both methods on calibrated visibilities obtained on the 
same target on 4 nights of very different seeing conditions. We note that the night-to-night varia- 
tions observed in this example are not believed to be due to the object itself (AB Aur) , but are 
instrumental. Variations in the object visibility are expected if (1) the source is asymmetric on 
the sky, or (2) it is composed of more than one component with changing relative fluxes. However, 
during one of the nights this source has been monitored for 5.2 hours (solid bullets in the plot), 
showing no sign of variation and implying a high degree of circular symmetry. Moreover, according 
to Voshchinnikov, Molster, & The (1996), AB Aur is an irregular photometric variable of very small 
amplitude. Am = ±0.02. In a representation of the source consisting of a central star surrounded 
by near-IR emitting material such a small change in the brightness of the central star would result 
in a change of the observed visibility < 1%, and therefore this effect is not likely to be the cause of 
the observed night-to-night variations. 
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Fig. 2. — Calibrated visibilities (Vtme) obtained using the visibility distribution model (top panel) 
and the mean of the measured distribution (bottom panel). Each symbol represents data taken 
on the same source (AB Aur) on different nights covering a large range of seeing conditions, and 
using identical detection parameters. For reference, the longest set of data (filled bullets) spans 5.2 
hours. 
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3. THE SAMPLE, LOG OF OBSERVATIONS, AND VISIBILITY AMPLITUDES 

To date, the most complete and comprehensive compilation of candidate HAEBE stars is the 
catalog of 287 objects by The, de Winter, &; Perez (1994). Our sample was selected from their 
Table 1, which contains 108 Ae and Be stars recognized as true members or good candidates of 
the HAEBE group, using as the selection criteria that the stars be brighter than the limiting 
magnitudes of lOTA's tip-tilt servo system (V = 11 — 12) and ncar-IR fringe detector (H=7, 
K'=6.2). Furthermore, in order to have as homogenous a sample as possible, and in order to test 
the accretion disk hypothesis, we have primarily selected our sources to belong to Group I in the 
classification of HSYK. The only exceptions are u Ori and MWC 166, which are Group III stars in 
that classification. 

The final list of targets observed is presented in Table 3.2. Columns 1 — 3 contain the star 
names and coordinates; columns 4 — 10 show the VRIJHK photometry compilation, mostly from 
HSYK; columns 11 — 15 contain stellar parameters, distance, visual extinction, spectral type and 
effective temperature, mostly from the Hipparcos data study of van den Ancker et al. (1998); 
and the last column comments on whether the star posseses a visual or spectroscopic companion 
and/or IR excess. Table 2 shows a log of the observations, including the calendar date, spectral 
band, calibrator used, the approximate hour angle (HA) range covered during the observation, a 
qualitative measure of the seeing at the time of the observations (l=very good, 4=bad) based on 
real-time visual inspection of the amount of atmospherically induced fringe piston and amplitude 
noise, and the baseline used (identified by its approximate length). 

3.1. Special Cases and Visibility Correction Factors 

3.1.1. Stars with Close Visual Companions 

Using the technique of adaptive optics in the near-IR, Corporon (1998) surveyed 68 HAEBE 
stars and identified 30 binary systems with separations greater than 0.1 arcsec. We consider here the 
effect of those "wide" companions on our data set. Indeed, if a star has a companion at an angular 
offset such that it is in the detector's field of view (~ 3" x3" for each pixel of the NICM0S3 camera), 
but whose interferogram has a zero OPD point that does not fall near that from the primary 
source, in the resulting scan a correction needs to be made to the fringe visibility measured on 
each interferogram that accounts for the incoherent flux added by the other star (exactly analogous 
to the incoherent flux added by thermal background, which needs to be subtracted before the 
fringe amplitude is measured). If we call V^'J.^g the calibrated visibility corresponding to the bright 
(j = ^4) or faint {j = B) fringe obtained by the procedure described in Section 2.2, then the 
corrected visibility is given by 



y'true = yLe/rj\ with: rA,B = ^ iOTAm/2.5) 



(4) 
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where Am = vriB — ruA is the magnitude difference between the binary components in the relevant 
spectral band. 

Our sample overlaps almost perfectly with that of Corporon (1998), and therefore we use 
his measured binary separations and flux ratios to apply the above correction in the appropriate 
cases. Table 3 summarizes the angular separations (9), magnitude differences and the spectral types 
assigned to each component, for the subset of our sample for which the adaptive optics companion 
is approximately within the field of view of our NICM0S3 camera. Note that in our data set, 
the interferogram corresponding to the faint companion is never detected, either because it is too 
faint or because it occurs at a value of the OPD which is outside of the scan length. However, the 
correction by the factors, last two columns in Table 3, still needs to be applied to the fringe 
visibility corresponding to the bright component. The only two stars in our sample which were not 
observed by Corporon are AB Aur and V594 Cas, and we will assume here that they are single 
stars. 



3.1.2. MWC 1 080 and HR8881 

The data for MWC 1080 shows that it is more resolved at H than at K', by an amount signif- 
icantly greater than expected from the higher resolving power of the shorter H-band wavelength. 

In searching for explanations for this unusual circumstance, it was discovered that according to 
Hammcrslcy ct al. (1994), the calibrator used for MWC 1080 (HR 8881), is misclassified in the 
Bright Star Catalogue^, and is really a M3 giant or K supergiant, and not a Kl dwarf. 

Assuming that HR 8881 is a normal giant, we may estimate its angular size and account for it 
in the calibration of the MWC 1080 data. Based on a compilation of angular diameters measured 
in the near-IR by the techniques of lunar occupations and long baseline inteferometry, van Belle 
(1999) has derived empirical relations between angular diameters and {V — K) and {B — K) colors 
for main sequence, giants and supergiants stars. Using the photometry and Ay = 0.22 estimate 
of Hammersley et al. (1994), we obtain {V - K) = 4.05 ± 0.02 and {B - K) = 5.33 ± 0.02. The 
empirical relations using those two colors predict uniform disk (UD) diameters 6 = 2.12 it 0.28 mas 
and 6 = 2.08 it 0.37 mas respectively. The weighted mean of the two predictions is ^ = 2.10 ± 0.22 
mas, which we adopt as our estimate of the angular size of HR 8881. 

The visibility data for MWC 1080 has therefore been calibrated using the visibility that is 
expected for the estimated UD diameter V^f = 2 ■ JiiirOs)/ {-kOs). For reference, for the baseline 
lengths corresponding to those obervations, this calibration factor is on average 0.95 for the H band 
observations and 0.97 for the K' band observations. 



^The Bright Star Catalogue, Fourth Revised Edition, Yale University Observatory, 1982 
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The unusual circumstance therefore remains that MWC 1080 is resolved by about 20% in H 
band observations and appears essentially unresolved in K' observations. This difference is difficult 
to explain with any simple model of the source, given that the near-IR excess is greater at K than at 
H, contributing 79% and 91% of the total flux at H and K respectively (see Section 4). One possible 
explanation would be that the calibrator, HR 8881, is not a UD, and that the interferometer resolves 
it at K' by a greater amount than is predicted by its photospheric extension. There are however 
no specific measurements in the literature to confirm or rule out this hypothesis, and a complete 
resolution to this question will have to await new observations at the IOTA of both MWC 1080 
and HR 8881. If the result is confirmed, the apparent decrease in size with increasing wavelength 
could be explained if MWC 1080 is an under-resolved core-halo object. In this model, the small 
core contains the central star as well as the source of near-IR excess and dominates the emission 
at K' wavelengths, while at the shorter H wavelengths scattering reveals the larger angular extent 
of the cold halo component. 



3.2. Visibility Data 

The final set of calibrated H and K'-band fringe visibility data is shown in Figure 3. The 
data are presented in two panels, one for each independent variable specifying the projection of the 
baseline vector in a plane (u, v) normal to the direction toward the source. The left panels show the 
data as a function of the modulus of the projected baseline: s = Vu'^ + v'^, expressed in millions of 
wavelengths at the center of the observation band. The left panels show the data as a function of 
the projected baseline position angle: $ = arctan (u/'u), measured from North toward East. The 
size of the data symbols equals our typical uncertainty of 5% in the calibrated visibilities. In the 
case of the stars known to be adaptive optics binaries, only the bright component (A) detected 
interferometrically appears in the figure, and explicit error bars are shown reflecting the error in 
correction by the flux factor described in Section 3.1.1. The error bars in the data for MWC 1080 
also include the contribution due to the uncertainty in the estimated size of HR 8881. 

The detailed interpretation of this data set is the subject of the next two sections. However, 
based on inspection, we may immediately make the following observations: (1) of the 15 stars 
observed, 11 (73% of the sample) are resolved by the interferometer; (2) one star, MWC 361- A, 
shows the clear sinusoid signature of a binary system; (3) none of the other stars shows significant 
variation of visibility with baseline position angle (the range covered by the observations varies from 
A$ = to 65°); and (4) of the two Group HI sources, to Ori appears unresolved, and MWC 166 
appears only slightly resolved, both consistent with that classification. 



Table 1. Target list 



Name 


RA(2000) 


DEC(2000) 


V 


R 


I 


J 


H 


K 


Ref. 


d(pc) 


A*v 


SpTy 


Log(Teff) 


Ref. 


Comments 


Other Names 


V594 Cas 


00 


43 


: 18.0 


+61 


: 54 : 


37.0 


10.38 


9.92C 


9.40C 


8.50 


6.90 


5.70 


(1) 


(650) 


2.14 


B8eq 


4.08 


(2) 


IRE 


MWC 419 


AB Aur 


04 


55 


: 45.8 


+30 


: 33 : 


05.5 


7.07 


6.82J 


6.69J 


6.10 


5.10 


4.40 


(1) 


144 


0.50 


AOVe 


3.98 


(2) 


IRE 


MWC 93/HD 31293 


T Ori 


05 


35 


: 50.0 


-05 


: 28 : 


42.0 


10.35 


10.02C 


9.65C 


8.45 


7.27 


6.22 


(1) 


460 


1.7 


B9e 


4.03 


(1) 


VB/SB/IRE 


MWC 763 


V38Q Ori 


05 


36 


: 25.3 


-06 


: 42 : 


57.1 


10.37 


9.95C 


9.43C 


8.21 


7.18 


6.09 


(1) 


(460) 


1.43 


Ale 


3.97 


(2) 


VB/SB/IRE 


MWC 765 


u} Ori 


05 


39 


: 11.1 


+04 


: 07 : 


17.2 


4.57 


4.57J 


4.67J 


4.72 


4.73 


4.71 


(1) 


(400) 


0.28 


B3nie 


4.23 


(2) 




MWC 117 


MWC 147 


06 


33 


: 05.1 


+10 


: 19 : 


20.0 


8.70 


8.40C 


8. IOC 


7.48 


6.67 


5.68 


(1) 


(800) 


1.61 


BlVe 


4.41 


(2) 


VB/IRE 


HD 259431 


MWC 166 


07 : 


04 


: 25.5 


-10 


: 27 : 


16.7 


7.01 


6.65C 


6.24C 


5.79 


5.53 


5.24 


(1) 


(1150) 


2.23 


BOIVe 


4.50 


(2) 


VB/SB 


HD 53367 


MWC 863 


16 : 


40 


: 17.8 


-23 


: 53 : 


44.1 


8.80 


8.41C 


7.97C 


7.06 


6.32 


5.51 


(1) 


150 


1.61 


AlVe 


3.97 


(2) 


VB/SB/IRE 


HD 150193 


MWC 275 


17 : 


56 


: 21.3 


-21 


: 57 : 


19.5 


6.83 


6.77C 


6.72C 


6.25 


5.57 


4.72 


(1) 


122 


0.25 


AlVe 


3.97 


(2) 


IRE 


HD 163296 


MWC 297 


18 : 


27 


: 39.5 


-03 


: 49 : 


52.0 


12.17 


9.45J 


8.39J 


6.01 


4.58 


3.08 


(1) 


450 


8.3 


09e 


4.53 


(1) 


IRE 




MWC 614 


19 : 


11 


: 11.1 


+15 


: 47 : 


16.6 


8.67 








6.7 


5.91 


(3) 


240 


1.27 


B9e 


4.02 


(2) 


IRE 


HD 179218 


V1295 Aql 


20 : 


03 


: 02.4 


+05 


: 44 : 


17.4 


7.84 








6.57 


5.78 


(3) 




0.19 


A2nie 


3.95 


(2) 


IRE 


MWC 325/HD 190073 


V1685 Cyg 


20 : 


20 


: 28.1 


+41 


: 21 : 


53.0 


10.54 






7.90 


6.91 


5.78 


(1) 


(980) 


3.16 


B2Ve 


4.34 


(2) 


VB/IRE 


MWC 340 


MWC 361 


21 : 


01 


: 36.8 


+68 


: 09 : 


48.1 


7.36 


6.88J 


6.35J 


6.10 


5.50 


4.70 


(1) 


430 


1.92 


B2IVe 


4.31 


(2) 


VB/SB/IRE 


HD 200775 


MWC 1080 


23 : 


17 


: 26.0 


+60 


: 50 : 


43.0 


11.68 


10.55J 


8.92J 


7.40 


6.00 


4.70 


(1) 


(2200) 


5.27 


BOe 


4.48 


(2) 


VB/SB/IRE 





References. — (1) Hillenbrand et al. 1992 (2) van den Ancker et al. 1998 (3) Malfait et al. 1998 

Note. — 

(a) J=Joiinson, C=CoLisins. 

(b) Distances in curved brackets are tlie values from tlie literature adopted in van der Ancker et al. 1998. 

(c) H band magnitude of MWC 614 (6.7±0.2) estimated from IOTA data. 

(d) VB— Visual Binary (Corporon 1998), SB— Spectroscopic Binary (from Corporon Kl Lagrange 1999), IRE— IR Excess (from Hillenbrand et al. 1992). 
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Table 2. Log of Herbig Ae/Be Stars Observations 



Date 


Target 


Band 


Calibrator 


HA 


Seeing^ 


Baseline 








(HR) 






(m) 


1997 Oct 13 


AB Aur 


H 


1343,2219 


[-1.2, +0.6] 


2 


21 


1997 Nov 15 


AB Aur 


TT 

H 


1626 


[—1.0, +3.0] 


1 


O CI 

38 


1998 Jan 19 


AB Aur 


K 


1626 


[+0.9, +1.1J 


2 


O CI 

38 


1998 Mar 02 


AB Aur 


K 


1626 


[+1.9, +2.1] 


1 


21 


1998 Mar 03 


uj On 


H 


2019 


[+0.3, +0.7J 


2 


21 






K 


2019 


[+0.9, +1.1] 


2 


21 


1998 Mar 07 


A T ~> A 

AB Aur 


H 


1626 


[+1.6, +1.9] 


2 


38 






K 


1626 


r 1 r\ 1 p' T 

[+2.2, +2.5] 


2 


38 


1998 Jun 12 


MWC 297 


H 


7149 


[+0.9, +1.2] 


3 


38 




MWC 361 


H 


7967 


[—0.4, +0.2] 


3 


38 


1998 Jun 13 


MWC 863 


TT 

H 


6153 


[—0.8, —0.5] 


2 


O CI 

38 




MWC 275 


H 


6704 


[—0.9, —0.1] 


2 


38 




MWC 361 


H 


7967 


[+0.4, +0.7] 


1 


38 


1998 Jun 14 


MWC 297 


TT 


7149 


r n o r\ r\\ 

[—0.8, —0.0] 


2 


38 


1998 Jun 15 


MWC 361 


H 


7967 


[—0.5, —0.0] 


1 


38 


1998 Jun 17 


MWC 863 


H 


6153 


[—0.1, —0.5] 


1 


38 




MWC 275 


H 


6704 


[—0.3, +1.4] 


1 


38 




MWC 361 


TT 

H 


7967 


[—0.5, +0.7] 


1 


38 


1998 Jun 18 


MWC 297 


TT 

H 


7066 


[—1.9, —1.6] 


1 


o o 

38 




MWC 361 


TT 


7967 


r r\ o rv oT 

[—0.8, —0.2] 


1 


38 






K 


7967 


[+0.0, +0.2] 


1 


38 




MWC 1080 


TT 


8881 


[—1.8, —1.4] 


1 


o o 

38 


1998 Jun 19 


MWC 361 


H 


7967 


[-2.8, +1.0] 


1 


38 






K' 


7967 


[-2.1, +0.4] 


1 


38 


1998 Jun 20 


MWC 361 


H 


7967 


[-4.9, +0.4] 


1 


38 






K' 


7967 


[-4.3, +0.6] 


1 


38 




MWC 1080 


H 


8881 


[-3.5,-2.9] 


1 


38 


1998 Sep 28 


MWC 361 


H 


7967 


[-0.7, +4.4] 


1 


38 






K' 


7967 


[+2.1, +3.7] 


1 


38 


1998 Sep 29 


V1295 Aql 


H 


7669 


[+1.0, +2.4] 


2 


38 


1998 Sep 30 


MWC 614 


H 


7449 


[+1.2, +2.4] 


1 


38 


1998 Oct 02 


MWC 614 


H 


7449 


[+1.3, +2.2] 


1 


38 
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Table 2 — Continued 



Date 


Target 


Band 


Calibrator 
(HR) 


HA 


Seeing^ 


Baseline 
(m) 




V1685 Cyg 


H 


7800 


[+2.4, +3.7] 


2 


38 


1998 Oct 31 


MWC 1080 


H 


8881 


[+0.9, +2.9] 


1 


38 






K' 


8881 


[+1.8, +3.8] 


1 


38 




T Ori 


H 


1986 


[-0.0, +0.4] 


1 


38 


1998 Nov 02 


AB Aur 


H 


1626 


[-2.3, +2.9] 


1 


38 


1998 Nov 03 


V594 Cas 


H 


237 


[-1.9, +0.8] 


1 


38 




V380 Ori 


H 


1697 


[-1.8,-0.4] 


1 


38 




MWC 147 


H 


2426 


[-0.2, +0.4] 


1 


38 






K' 


2426 


[+0.8, +1.3] 


1 


38 


1998 Nov 07 


u) Ori 


H 


2019 


[-1.4, +2.0] 


3 


38 


1998 Nov 10 


MWC 1080 


K' 


8881 


[+0.0, +3.0] 


2 


38 




MWC 147 


H 


2426 


[-2.5,-0.0] 


2 


38 




MWC 166 


H 


2723 


[+0.1, +1.5] 


3 


38 


1999 Mar 01 


V380 Ori 


H 


1986 


1.9 


2 


21 




MWC 863 


H 


6153 


[-0.8,-0.4] 


2 


21 


1999 Mar 02 


MWC 863 


H 


2723 


[-2.3,-0.9] 


2 


21 






K' 


2723 




2 


21 


1999 Mar 05 


MWC 166 


H 


2723 


[+1.7, +2.2] 


2 


38 



^The qualitative measure of seeing goes from "good" 



= 1 to "bad" = 3 
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Table 3. HAEBE stars with adaptive optics companions 



Name 


e" 


AV 


AR 


AI 


AJ 


AH 


AK 


SpTy 


SpTy 


rA{H) 


rA{K') 


















A 


B 






V380 Ori 


0.14 




0.10 




0.11 


0.94 


1.36 


B7 


GO 


0.704 ± 0.04 


0.778 ± 0.03 


MWC 147 


3.1 


7.22 


6.82 


6.02 


4.81 


5.05 


5.67 


B2 


K5 


0.990 ± 0.002 


0.995 ± 0.001 


MWC 166 


0.64 


1.43 


1.61 


1.7 


1.70 


1.40 


1.77 


BO 


BO 


0.784 ± 0.03 


0.836 ± 0.02 


MWC 863 


1.07 


3.57 




2.57 


2.07 


2.24 


2.72 


A2 


K4 


0.887 ± 0.02 


0.924 ± 0.01 


V1685 Cyg 


0.72 












5.50 


B2 






0.994 ± 0.001 


MWC 361 


2.55 


5.40 


5.04 


4.52 




4.12 


4.64 


B2 


G5 


0.978 ± 0.004 


0.986 ± 0.002 


MWC 1080 


0.77 


2.34 


2.44 


2.61 




2.84 


2.28 


B2 


B2 


0.932 ± 0.01 


0.891 ± 0.02 



References. — Angular separation, photometric and spectral type data from Corporon (1998). 

Note. — (1) Completeness limits: 0.1" <0< 7"; Am < 6 - 7; (2) Uncertainties: ae = 0.05"; aAm < 0.2. 
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Fig. 3. — H and K'-band calibrated visibility data. For each source, the data are plotted as a 
function of the projected baseline length (s, left panel) and position angle (<J>, right panel). The 
different symbols correspond to observations on different nights, and retain their meaning between 
the right and left panels for a given source. Observations on the 21 m baseline are represented 
with filled symbols, and observations at the 38 m baseline with open symbols. All subsequent plots 
involving this visibility data retain the same symbols. 
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4. EXCESS NEAR-INFRARED FLUX 

Most of the HAEBE stars in this paper are known to have a strong near-IR excess, and are 
resolved by the interferometer. Therefore, in order to interpret the visibihty amplitude data we 
consider a two component model consisting of the central star, plus a component which is the source 
of the near-IR excess. However, before this model can be used to interpret the interferometer data, 
the relative contributions to the total H and K fluxes from these two components must be quantifled. 
The notation used hereafter is as follows: "T", 'V and superscripts denote total, stellar and 
excess fluxes respectively. Primed fluxes are apparent, while un-primed fluxes are intrinsic (i.e. 
corrected for extinction). 

4.1. Stellar Flux 

In order to estimate the stellar near-IR fluxes, we use the photometry data, visual stellar 
extinctions and stellar eflFective temperatures from the literature tabulated in Table 3.2. We assume 
that the short wavelength (VRI) flux arises entirely from the central star, and approximate the 
stellar spectrum as a blackbody at the adopted effective temperature. The solid angle subtended 
by the stellar disk of radius i?* at distance d, = TxR^/d"^, is the parameter adjusted in order to 
provide a best fit to the de-reddened V, R and I fluxes, according to: -Fy^j = BvRi{Teff) ■ fi^, 
where Bi,(T) denotes the Planck function. The literature flux data has been de-reddened using 
the extinction law [Ax/Av) of Steenman & The (1991) and a ratio of total to selective extinction 
R = Ay/Eb-v = 3.10. Once the solid angle factor is thus determined, the stellar flux may be 
calculated at any frequency as F* = B^{Teff) ■ il*, and used to quantify the flux in excess of the 
photospheric value. 

The literature flux data is quoted as having uncertainties of 0.01 — 0.02 mag in the VRI bands. 
In addition, the spectral classification is uncertain by 1 sub-class typically, which results in about 
0.2 mag uncertainty in the derived values of the visual extinction and about 11% uncertainty in the 
effective temperatures. The resulting relative uncertainties in the predicted near-IR stellar fluxes 
are in the range 17% - 50% (median 30%). 

In the case of those sources which were identified as adaptive optics binaries by Corporon, we 
calculate the SED of each component by combining the magnitude differences of Table 3 with the 
total photometry of Table 3.2. We determine effective temperatures for each component using the 
spectral types assigned by Corporon and the temperature scale of Cohen & Kuhi (1979). These 
effective temperatures are used to calculate the blackbody model which represents the photospheric 
contribution of each star to the total fiux, as described above for the single stars. 
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4.2. Excess Flux 

Using this model to represent the stellar photosphere, we may deduce the apparent excess fluxes 
due to CS emission: Fj^^^ = F'J^^—Fj^^ ■ 1{)''^1jhkI'^-^ ^ where the apparent total fluxes and visual 
stellar extinctions are the observed quantities from the literature. Estimating the intrinsic (i.e. de- 
reddened) excess fluxes is also of interest, since they directly relate to the physical properties of 
the source. This however, requires knowledge of the amount of attenuation suffered by the excess 
radiation between the source and the observer, while the only information available is the total 
extinction in the line of sight toward the star. Therefore, we adopt here the simplifying assumptions 
that (1) all of the stellar extinction is due to cold foreground material, and (2) it is the same for 
aU lines of sight that intercept the extended component, so that Ffjjj^ = 10+^ 

Figure 4 shows the de-reddened SED data and best fit stellar blackbody models for all the 
target stars. Note that in most cases the individual components in the adaptive optics binaries also 
have significant IR excess. 

The results for the derived H and K stellar fluxes are summarized in Table 4. In addition to 
the quantities already deflned, the fourth column lists the stellar angular diameters that are implied 
by the solid angle calculated, 6* = 2-\/u/tv, and it can be seen that the values obtained are all 
small compared to the resolution of the interferometer. Therefore, in our models of the brightness 
distribution, the central stars will be treated as point sources. 

Tabic 5 summarizes the results obtained for the excess fluxes. The errors are dominated 
by the uncertainty in the determination on the underlying stellar spectrum, as described in the 
previous section, and include an uncertainty of 0.01 — 0.1 mag from the total near-IR magnitudes 
in the literature, as well as a small contribution due to the uncertainty in the extinction at near-IR 
wavelengths derived from the visual extinctions. The last column in Table 5 contains the color 
temperature (Tc) of the excess emission, calculated from the intrinsic H and K excesses such that 
F§/F§ = Bh{Tc)/Bk{Tc). 
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Fig. 4. — Total SED data and photospheric models. Each panel shows the de-reddencd total 
SED data of a target (circles) and the stellar blackbody model (dashed lines) corresponding to the 
underlying star. In the case of the adaptive optics binaries (Corporon 1998), the SED data for each 
component is also shown (open squares and triangles) and the solid line represents the sum of the 
two model photospheres. 
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Table 4. Near-IR stellar fluxes 



Name 


log(Teff) 




6** 








Ah 


Ak 














(10-2° sr) 


(mas) 


(Jy) 


(Jy) 






(Jy) 


(Jy) 


V594 Gas 


4.08 


6.3 


0.06 


0.5 


± 


0.1 


0.3 


± 


0.1 


0.3 


0.2 


0.4 


± 


0.1 


0.3 


± 


0.1 


AB Aur 


3.98 


50.3 


0.2 


3.0 


± 


1.1 


1.9 


± 


0.7 


0.07 


0.04 


2.8 


± 


1.0 


1.8 


± 


0.7 


T Ori 


4.03 


5.4 


0.05 


0.4 


± 


0.1 


0.2 


± 


0.1 


0.2 


0.1 


0.3 


± 


0.1 


0.2 


± 


0.1 


V380 Ori-A 


4.09 


2.2 


0.03 


0.2 


± 


0.05 


0.1 


± 


0.03 


0.2 


0.1 


0.1 


± 


0.04 


0.1 


± 


0.03 


to Ori 


4.23 


98.0 


0.2 


12.9 


± 


3.8 


7.8 


± 


2.4 


0.04 


0.02 


12.4 


± 


3.7 


7.6 


± 


2.3 


MWC 147-A 


4.31 


6.0 


0.06 


1.0 


± 


0.2 


0.6 


± 


0.1 


0.2 


0.1 


0.8 


± 


0.2 


0.5 


± 


0.1 


MWC 166-A 


4.49 


22.0 


0.05 


6.0 


± 


1.3 


3.5 


± 


0.8 


0.3 


0.2 


4.2 


± 


0.9 


2.9 


± 


0.6 


MWC 863-A 


3.96 


30.9 


0.1 


1.7 


± 


0.6 


1.1 


± 


0.4 


0.2 


0.1 


1.3 


± 


0.5 


1.0 


± 


0.3 


MWC 275 


3.97 


51.2 


0.2 


2.9 


± 


1.1 


1.9 


± 


0.8 


0.03 


0.02 


2.8 


± 


1.1 


1.8 


± 


0.7 


MWC 297 


4.53 


60.7 


0.2 


18.3 


± 


3.4 


10.6 


± 


2.0 


1.2 


0.7 


6.3 




1.2 


5.8 


± 


1.1 


MWC 614 


4.02 


19.6 


0.1 


1.3 


± 


0.4 


0.8 


± 


0.3 


0.2 


0.1 


1.1 


± 


0.4 


0.8 


± 


0.2 


V1295 Aql 


3.95 


24.9 


0.1 


1.3 


± 


0.5 


0.9 


± 


0.3 


0.03 


0.01 


1.3 


± 


0.4 


0.8 


± 


0.3 


V1685 Cyg 


4.34 


4.1 


0.05 


0.7 


± 


0.2 


0.4 


± 


0.1 


0.4 


0.2 


0.5 


± 


0.1 


0.4 


± 


0.1 


MWC 361- A 


4.31 


29.5 


0.1 


4.9 


± 


1.2 


2.9 


± 


0.7 


0.3 


0.1 


3.6 


± 


0.9 


2.5 


± 


0.6 


MWC 1080-A 


4.31 


10.4 


0.07 


1.7 


± 


0.4 


1.0 


± 


0.2 


0.7 


0.4 


0.8 


± 


0.2 


0.7 


± 


0.2 



Table 5. Near-IR excess fluxes 



Name Fif F';[ F'i F'i {F'^/F'^)h {F'^/F'^)k F§ F§ n 

(Jy) (Jy) (Jy) (Jy) (%) (%) (Jy) (Jy) (k) 



V594 Cas 


1.7 


± 


0.1 


3.3 


± 


0.2 


1.3 


± 


0.1 


3.0 


± 


0.2 


76.2 


± 


7.0 


91.5 


± 


2.7 


1.7 


± 


0.2 


3.5 


± 


0.2 


1382 


± 


117 


AB Aur 


8.9 


± 


0.4 


10.8 


± 


0.5 


6.1 


± 


1.1 


8.9 


± 


0.9 


68.7 


± 


11.5 


82.9 


± 


6.5 


6.5 


± 


1.2 


9.3 


± 


0.9 


1814 


± 


331 


T Ori 


1.2 


± 


0.01 


2.0 


± 


0.1 


0.9 


± 


0.1 


1.8 


± 


0.1 


74.8 


± 


8.6 


89.5 


± 


3.7 


1.1 


± 


0.1 


2.0 


± 


0.1 


1503 


± 


110 


V380 Ori-A 


0.9 


± 


0.05 


1.8 


± 


0.1 


0.8 


± 


0.1 


1.7 


± 


0.1 


84.3 


± 


4.3 


94.2 


± 


1.6 


1.0 


± 


0.1 


1.9 


± 


0.1 


1457 


± 


113 


w Ori 


12.6 


± 


0.6 


8.1 


± 


0.4 


0.1 


± 


3.7 


0.5 


± 


2.4 


0.9 


± 


29.8 


5.9 


± 


29.1 


0.1 


± 


3.8 


0.5 


± 


2.4 








MWC 147-A 


2.1 


± 


0.1 


3.3 


± 


0.2 


1.3 


± 


0.2 


2.8 


± 


0.2 


62.8 


± 


9.2 


84.3 


± 


3.9 


1.7 


± 


0.3 


3.2 


± 


0.2 


1467 


± 


189 


MWC 166-A 


4.7 


± 


0.2 


4.2 


± 


0.2 


0.5 


± 


0.9 


1.3 


± 


0.7 


11.7 


± 


20.1 


31.2 


± 


15.8 


0.8 


± 


1.4 


1.6 


± 


0.8 


1408 


± 


1690 


MWC 863-A 


2.6 


± 


0.1 


3.6 


± 


0.2 


1.2 


± 


0.5 


2.6 


± 


0.4 


47.8 


± 


18.2 


73.1 


± 


9.8 


1.6 


± 


0.6 


3.0 


± 


0.4 


1471 


± 


405 


MWC 275 


5.8 


± 


0.3 


8.0 


± 


0.4 


3.0 


± 


1.1 


6.2 


± 


0.8 


51.1 


± 


19.0 


77.0 


± 


9.3 


3.1 


± 


1.2 


6.3 


± 


0.9 


1393 


± 


375 


MWC 297 


14.4 


± 


0.7 


36.3 


± 


1.8 


8.2 


± 


1.4 


30.6 


± 


2.1 


56.6 


± 


8.5 


84.1 


± 


3.1 


23.8 


± 


4.2 


56.3 


± 


4.1 


1268 


± 


143 


MWC 614 


2.0 


± 


0.1 


2.7 


± 


0.1 


0.9 


± 


0.4 


1.9 


± 


0.3 


44.7 


± 


17.8 


71.4 


± 


9.6 


1.1 


± 


0.4 


2.1 


± 


0.3 


1453 


± 


387 


V1295 Aql 


2.3 


± 


0.1 


3.0 


± 


0.1 


1.0 


± 


0.5 


2.2 


± 


0.3 


44.0 


± 


19.8 


72.0 


± 


10.3 


1.0 


± 


0.5 


2.2 


± 


0.3 


1325 


± 


424 


V1685 Cyg 


1.7 


± 


0.1 


3.0 


± 


0.1 


1.2 


± 


0.1 


2.7 


± 


0.2 


70.6 


± 


7.2 


88.4 


± 


2.9 


1.8 


± 


0.2 


3.4 


± 


0.2 


1463 


± 


126 


MWC 361-A 


6.0 


± 


0.3 


8.1 


± 


0.4 


2.4 


± 


0.9 


5.6 


± 


0.7 


40.5 


± 


14.8 


69.4 


± 


7.7 


3.3 


± 


1.3 


6.6 


± 


0.9 


1408 


± 


387 


MWC 1080-A 


3.6 


± 


0.2 


7.3 


± 


0.4 


2.9 


± 


0.3 


6.6 


± 


0.4 


79.2 


± 


5.2 


91.0 


± 


2.3 


6.6 


± 


0.6 


10.4 


± 


0.7 


1673 


± 


151 
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5. MODELS OF THE BRIGHTNESS DISTRIBUTION 

The small sizes of the central stars derived in Section 4 imply that the interferometer resolves 
the CS material source of the excess emission. In this section we make use of the interferometer data 
to characterize the geometrical distribution of this CS material. This information, together with 
the intrinsic excess fluxes derived in Section 4, permits the brightness distribution to be estimated. 

5.1. A Source with no IR Excess: uj Ori 

The visibility data of Figure 5 shows that u) Ori is unresolved by the interferometer. More- 
over, it was shown in Section 4 that there is no significant IR excess associated with this star. 
Consequently, we will place a size limit to the star itself by considering a single uniform disk (UD) 
component. 

Using the visibility data obtained at the longest 38 m baseline, which provides the highest 
angular resolution, the mean fringe spacing projected on the sky is s = 21. IMA and the mean 
visibility is y = 100.5 ± 2.0%, where the error is the error in the mean of the 5 long baseline 
measurements. Therefore the minimum visibility allowed by the data is Vmin = 96.5% and by 
solving in the equation for the UD visibility (recall Section 3.1.2) it is found that the upper limit 
to the UD angular diameter is 9 < 1.3 mas, at the 2a level. The data and visibility function 
corresponding to this upper limit are shown in Figure 5. 

We note that among our sample, this is the source which most clearly shows the absence of IR 
excess above photosphcric levels (sec Figure 4), and as such was classified as Group III by HSVK. 
Therefore, observations of this source are valuable as a control experiment, and our result is indeed 
consistent with it being a naked photosphere. 




10 20 30 

Spacing (MA) 

Fig. 5. — uj Ori visibility data (H-band) and model curve corresponding to the upper limit UD 
diameter {9 < 1.3 mas). 
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5.2. Models of the Source of IR Excess 

For the remainder of our sample, we consider four plausible models of the source brightness. 
Each model has two components, with the stellar emission arising in a central point source and the 
excess near-IR emission contributed by either: (1) a Gaussian brightness distribution, (2) a uni- 
formly bright ring, (3) an infrared companion or (4) a "classical'" accretion disk with a temperature 
law T(r) oc r"'^/*^. The Gaussian and ring models are intended to provide a characteristic size for 
the source of near-IR excess, and consistent with the observation that the visibilities remain con- 
stant with projected baseline orientation, we consider the simplest case that they appear circularly 
symmetric on the sky. 

In general, the source specific intensity (i.e. brightness) is given by 



Iu{^,v) = F:- 6(0,0) +I^{^,ri) 



(5) 



and the corresponding normalized visibility amplitude is calculated by application of the Van 
Cittert-Zernike theorem (see e.g. Goodman 1985) 



Vu{u,v) 



FT{/,(e,77)} 



p* A. pE 



F^ + FT{I^{(,V)} 



F* + 



(6) 



where in the above equations, (^, r/) are angular offsets with respect to the source center along the 
projected baseline coordinates [u, v) measured in wavelength units (see e.g. Thompson, Moran, &; 
Swenson 1986), and FT denotes the Fourier transform^. 

In addition to u Ori, we find 3 more sources which appear unresolved, namely V594 Gas, T Ori 
and MWG 147. However, since in those cases there is significant near-IR excess (see Figure 4), we 
will place appropriate upper limits to the extended component within the context of the models 
considered. The source MWG 166 appears to have low visibility {V = 88.2 ± 0.9 %). However, it 
was found in the previous chapter to have a marginally significant amount of near-IR excess flux, 
with 88 lb 20 % of the total H fiux arising in the central star. The data is therefore consistent 
with this source having an extended component which is completely resolved by the interferometer. 
Given the large error in the derived excess flux however, we cannot rule out the possibility that 
there is zero excess, so that all the emission comes from the star alone, or that the interferometer 
only partially resolves an extended component which contributes a larger fraction (up to 32%) of 
the total H flux. In this case then, lower limits will be given for the size of the excess region under 
the models mentioned above. 



*Note that in the above and following visibility equations, and in order to simplify the notation, we use intrinsic 
rather than apparent fluxes, due to the fact that under our assumption that the central star and excess source suffer 
the same extinction the apparent and intrinsic fluxes relative to the total are the same. The sizes derived, however, 
do not depend on this assumption. 



5.2.1. Gaussian Intensity 



We consider a circularly symmetric brightness distribution given by 



I^{^,V) = 4°-e ^ (7) 
so that 

F* + ■ e 41n2'' 



(8) 



where is the peak intensity and Oq is the angular FWHM in both sky coordinates and rj). 

Once the Gaussian FWHM is found by fitting Equation 8 to the visibility data, the peak 
intensity is obtained from the constraint that the flux in this component matches the measured 
excess flux 

F^ = j j Il-e d^dv = I^,-nG (9) 

where Q.g = 7r^6'Q/41n2, is the effective solid angle. In practice, we find it most informative 
to express the peak intensity in terms of the equivalent peak brightness temperature, defined by 
7° = B,y{T^). The brightness temperature can then be readily compared with the color temperature 
derived from the H and K excess fiuxes in order to determine whether this model is consistent with 
blackbody emission at the color temperature of the H and K excess. 

Table 6 summarizes the results obtained. In addition to the best fitting Gaussian FWHM, the 

table includes the RMS of the fit evaluated as ^xLn/(^ - 1) with = Eili(^/"*" - Vi^"'^^)^, 
the number of individual observations for each star (N), the linear diameter (D) corresponding to 
the measured angular diameter calculated using the distances of Table 3.2, and the peak brightness 
temperature (Tg) calculated as described above. For the two stars which were observed in both the 
H and K' bands and are resolved (AB Aur and MWC 863- A), it is found that the sizes in this model 
are, within the errors, the same, and the table shows in those cases the size that is derived from 
a simultaneous fit to both data sets. For the other stars, for which an independent measurement 
of the size in the K'-band is not available, the peak brightness temperature at K' is calculated 
using the size derived from fitting the H-band data. In all cases (as well as in the following tables 
corresponding the other models under consideration), the errors in the angular sizes refiect the 
uncertainty in the relative flux in each component, given in Table 5, rather than the much smaller 
formal error (of order 0.1 mas) associated with fitting a given set of visibility data. Figure 6 shows 
the data and models (solid lines) computed with the parameters of Table 6. 

We note that for AB Aur, a slightly better agreement is obtained by fitting the H and K'-band 
data sets separately. In that case we derive: 6q = 4.1 it 0.5 mas and 6q' = 4.6 it 0.2 mas, with 
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fit RMS of 3.4% and 5.9% respectively, to be compared to 3.5% and 9.6% obtained by fitting the 
model of Table 6 to the H and K' data separately. However, although formally a better fit, the 
different H and K'-band sizes are within their errors. In Millan-Gabct et al. (1999b), based on a 
smaller data set and using a different data reduction method, the opposite result was obtained, that 
the characteristic size for AB Aur appeared somewhat smaller in the K' band than in the H band; 
although the difference was again judged to be most likely insignificant. This important question 
will need to be settled through improved accuracy of our visibility calibration. 

By comparing the peak brightness temperatures of Table 6 with the color temperatures of 
Table 5, it is apparent that for most objects, the brightness temperature is lower than the color 
temperature. One possible physical scenario for the Gaussian model is as an approximation to the 
brightness resulting from a circumstellar environment in which the dust particles are distributed in 
an envelope around the central star. A possible interpretation then is that the near-IR emission is 
optically thin, in which case the observed color temperature depends on the physical temperature 
of the dust grains, their emission properties and the optical depth of the CS cloud. We now 
further explore the consequences of this observation under the simplifiying assumption that a single 
temperature characterizes the dust envelope. If the near-IR frequency dependence of the optical 
depth in the circumstellar environment of HAEBE stars follows the properties characteristic of 
interstellar dust [Ay ~ A~^'^, due to dust particles smaller than the wavelength of light), we may 
derive the physical temperature (T) of the grains and the peak optical depth by solving in 



where the ratio of H to K-band extinction has again been taken from Steenman 8z The (1991) with 
Rv = 3.1. 

On the other hand, from studies of the stellar extinction toward HAEBE stars, some authors 
find evidence for anomalous extinction laws {Ry — 5) indicative of the presence of a significant 
population of grains with sizes greater (radius a ~ 0.2 /xm) than the average in the interstellar 
medium (Djie & The 1978; Pezzuto, Strafella, & Lorenzetti 1997). In the limit that the dust grains 
become larger than the wavelength of light, the cross-section for absorption and scattering becomes 
independent of wavelength and results in neutral extinction (as in interplanetary dust grains, of 
a ~ 1 — 100 iJ,m). Therefore, we also compute the peak optical depth that results from making the 
assumption that the opacity is wavelength independent 



Case A: I Bk{T^) 

T-O / 



Ah/Ak = 1.75 



(10) 




(11) 



The results obtained under these two assumptions are summarized in Table 7. The "Case A" 
assumption results in an estimate of the physical temperature (T) of the dust grains and H-band 
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Table 6. Gaussian model parameters. 



Name 




Fit RMS 


N 


D 


H-band Tg 


K'-band T% 




(mas) 


(%) 




(AU) 


(K) 


(K) 


V594 Gas 


< 0.8 




6 


< 0.5 


> 1938 


> 2212^ 


AB Aur 


4.1 ± 0.5 


4.6 


47 


0.6 ± 0.1 


1359 ± 65 


1258 ± 63 


T Ori 


< 1.4 




2 


< 0.6 


> 1443 


> 1424*^ 


V380 Ori-A 


2.4 ± 0.1 


2.0 


8 


1.1 ± 0.05 


1208 ± 22 


1144 ± 20^ 


MWC 147-A 


< 1.0 




9 


< 0.8 


> 1765 


> 1884 


MWG 166-A 


> 4.5 




8 


> 5.2 


< 1003 


< 916^ 


MWG 863-A 


7.6 ± 1.7 


6.7 


18 


1.1 ± 0.3 


963 ± 44 


865 ± 21 


MWG 275 




2.8 


9 


6+°-2 


1175 ± 104 


1124 ± 104*^ 


MWG 297 


5.81SI 


2.0 


10 


^•"-•-0.2 


1500 ± 77 


1594 ± 92^^ 


MWG 614 


7.6_2.5 


2.0 


13 


1-8-0.6 


923 ± 74 


826 ± 70^ 


V1295 Aql 


5.4_i.4 


4.3 


4 


a 


991 ± 78 


910 ± 68^ 


V1685 Gyg 


2.6 ± 0.1 


2.8 


3 


2.5 ± 0.1 


1284 ± 30 


1235 ± 30*^ 


MWG 1080-A 


2.7 ± 0.2 


6.4 


20 


5.9 ± 0.2 


1566 ± 33 


1535 ± 35 



'^No distance information available 

''Galculated using the size derived from H-band visibility data 
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peak optical depth, as well as the visual extinction of the star due to the CS cloud that would 
be implied, calculated according to Ay = ^///0.14 = (t///2) • (2.5/0.141n 10). For the "Case B" 
assumption, the physical temperature equals the color temperature of Table 5, and we calculate 
the frequency independent peak optical depth {j^) and the corresponding visual extinction, Ay = 
(r/2) • (2.5/ In 10). Finally, for comparison, the table reproduces in the last column the measured 
value of the visual extinction toward the star obtained from the literature (see Table 3.2). 

Prom the results in the last two tables, we extract the following conclusions. First, the observed 
fluxes and our measured sizes are consistent with emission from dust envelopes having finite optical 
thickness (i.e. Tg < Tc), in 10 out of 13 cases, even at the position of the peak intensity. We note 
that under the other extreme assumption concerning the extinction of excess emission, i.e. that 
it is zero and this material causes all of the stellar extinction, the intrinsic excess fluxes would be 
given by columns 4 and 5 of Table 5, which are lower and would lead to even smaller values of 
the peak brightness and peak optical depths inferred. The unphysical result Tg > Tc is found for 
one of the resolved sources (MWC 297) and for two of the three unresolved sources (V594 Cas and 
MWC 147), inconsistently with the model assumptions, and is most likely a consequence of the 
inadequacy of the assumption that there is a single temperature throughout the envelope rather 
than a radial temperature gradient. Second, the assumption that the extinction properties in the 
envelope are interstellar (case A) leads to relatively large values of the visual extinction (median 
Ay = 5.8) in the line of sight toward the central stars. Under this assumption then, a special 
geometry that allows a clear line of sight to the star itself would need to be invoked, as has been 
proposed by Bcrrilli et al. (1992) and Hartmann, Kenyon, & Calvet (1993) in the context of SED 
modelling. Third, under the assumption that the extinction in the envelope is neutral (case B), the 
visual extinctions to the stars (median Ay = 0.2) are somewhat lower than the measured values. 
However, the values derived are lower limits, and therefore consistent with the observed ones, since 
they do not include an interstellar contribution, or a contribution from the larger particles that are 
likely to exist in a realistic distribution of sizes. 
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Tablc 7. Gaussian model temperature and peak optical depth. 



Case A Case B Measured 



Name 


T 

(K) 




Av 


{T = Tc) 


Av 


Av 


AB Aur 


1423 ± 61 


1.4 ± 0.1 


5.3 


0.2 ± 0.1 


0.1 


0.5 


V380 Ori-A 


1238 ± 23 


1.8 ± 0.1 


5.8 


0.3 ± 0.05 


0.2 


1.43 


MWC 863-A 


1104 ± 110 


0.4 ± 0.2 


1.5 


0.04 ± 0.02 


0.02 


1.61 


MWC 275 


1195 ± 100 


2.1 ± 0.3 


8.2 


0.4 ± 0.2 


0.2 


0.25 


MWC 297 


a 






a 




8.3 


MWC 614 


1074 ± 58 


0.3 ± 0.1 


1.2 


0.03 ± 0.01 


0.02 


1.27 


V1295 Aql 


1072 ± 86 


0.7 ± 0.1 


2.7 


0.1 ± 0.05 


0.06 


0.19 


V1685 Cyg 


1299 ± 30 


2.5 ± 0.1 


9.7 


0.6 ± 0.1 


0.3 


3.16 


MWC 1080-A 


1570 ± 65 


4.1 ± 0.1 


15.9 


1.2 ± 0.1 


0.6 


5.27 



^No solution possible with > 

Note. — The optical depths and temperatures in this table are rounded from 
the solutions to equations 10 and 11, the resulting numerical loss of accuracy is 
however within the qouted error bars, which result primarily from the uncertainties 
in the excess fluxes. 
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5.2.2. Uniform Ring 

In the previous section it was shown that under the Gaussian model, the effective sizes measured 
by the interferometer require that the emission be optically thin over most of the area of the source, 
in order to match the observed fluxes. An alternative solution is to consider optically thick emission 
at the observed color temperature, but reduce the emitting surface area by concentrating it in a 
narrow ring around the star. 

For a uniform ring of intensity Ii, = Bi,{Tc), and inner and outer angular diameters Oi and 62 
respectively, we have 



Vu{s) 



B,{Tc) if \ei<V^^+W <h02 



otherwise 



+ 



2s 



2s 



(12) 
(13) 



where for each value of the inner diameter, the outer diameter is constrained by the requirement 
that the flux from the ring matches the excess flux: 62 = V^fWT^ with n = F^/B^{Tc). The 
only parameter to be fit to the visibility data is then the ring inner diameter. 

The results are shown in Table 8, where in addition to the quantities already defined, the last 
two columns contain the linear ring inner and outer diameters. Figure 6 shows the data and models 
(dot-short dash lines) computed with the parameters of Table 8. We note that the source V594 Cas 
(which appears unresolved) can not be represented by this model, due to the fact that for any inner 
ring diameter, the outer diameter required to match the excess flux results in a source that would 
appear resolved. Similarly, the large excess flux for MWC 297 implies that any solution for the ring 
model has a large size, much larger than the characteristic size given by the Gaussian model. The 
data can formally be well fit only at the peak of the second lobe of the ring visibility function, and 
we do not consider this a likely interpretation of the data for this source. 
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Table 8. Ring model parameters. 



Name 


(mas) 


FitRMS 
(%) 


N 


(mas) 


^1 
(AU) 


(AU) 


V594 Cas 


b 














AB Aur 


4.2 ± 0.4 


5, 


.2 


47 


4.7 ± 0.4 


0.6 ± 0.05 


0.7 ± 0.05 


T Ori 


< 1.4 






2 


< 2.05 


< 0.6 


< 0.9 


V380 Ori-A 


2.5 ± 0.1 


2, 


.0 


8 


3.0 ± 0.1 


1.15 ± 0.50 


1.4 ± 0.05 


MWC 147-A 


< 0.4 






12 


< 2.0 


< 0.3 


< 1.6 


MWC 166-A 


> 9.5 


2, 


,5 


8 


> 9.6 


> 10.9 


> 11.0 


MWC 863-A 


8.21?:^ 


6. 


.2 


18 




1 9+0.4 
-•■•^-0.2 


1 O+0.4 
^•"J-0.2 


MWC 275 




3. 


.1 


9 




6+°-2 
^•"-0.1 


'-'•'-0.1 


MWC 297^= 


21.9+1-2 


4. 


.1 


10 


25.0+10 


99+0.5 


11.25+0-45 


MWC 614 


6.7 ± 1.5 


2, 


.3 


13 


6.9 ± 1.4 


1.6 ± 0.4 


1.7 ± 0.3 


V1295 Aql 


5.4_i.2 


4. 


.4 


4 


5.8_i.i 


c 


c 


V1685 Cyg 


2.7 ± 0.15 


2, 


.8 


3 


3.4 ± 0.15 


2.55 ± 0.15 


3.3 ± 0.15 


MWC 1080-A 


2.3 ± 0.3 


9, 


.0 


20 


3.6 ± 0.3 


5.1 ± 0.7 


7.9 ± 0.7 



^Derived from 9i and flux 

'^No solution can match flux and visibility data simultaneously, see text 
'^The only formal solution requires Oq, see text 

^No distance information is available 



-32- 




Fig. 6. — Visibility data and Gaussian and ring models. Solid lines: Point source plus Gaussian 
intensity model. Dot-short dashed lines: Point source plus uniform ring model. The models are 
calculated for the best fit parameters of Tables 6 and 8. 
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5.2.3. Accretion Disk 



Given that the stars in our sample which the interferometer resolves are all Group I sources 
in the classification of HSVK, we consider in this section whether such an accretion disk model is 
consistent with the interferometer data. 

We consider a flat, blackbody disk heated by a central star of effective temperature Tg^p and 
radius B^. Energy balance between stellar light absorbed at distance r by one side of the disk and 
subsequent re-radiation at the equilibrium temperature, results in a temperature law in the disk, 
for i?^/r ^ 1 given by (Priedjung 1985) 



1X1/4 .^^,3/4 



TreAr)=T.-[jj -[^-fj (14) 

The contribution to the disk temperature resulting from viscous heating in a steady state, 
optically thick, self-gravitating disk, also for -R*/r ^ 1, may be shown to be (Lynden-Bell & 
Pringle 1974) 



r.„M=(^^^y'*..-3/^ (15, 

where Mace is the accretion rate, a is the Stefan-Boltzman constant, and the other symbols have 
their usual meanings. 

Considering the reprocessing and accretion processes simultaneously results in a radial tem- 
perature law given by T{r) = \T^^p{r) + r^cc('')] convenient units, we have 



r(r) 



2.55 X 



1„-S(^)V. + 5.26X10'«(^) 



Ma 



10-5 



1/4 



-3/4 



' AU 



TiAU ■ r 



-3/4 
AU 



(16) 
(17) 



Therefore, in this model, if the properties of the central star are known, the disk temperature 
by stellar heating is completely specified; and the heating due to viscous accretion depends only on 
the mass accretion rate. In this paper, when fitting the SED and visibility data, we treat Tiau as 
the parameter which sets the disk temperature scale, independent of the details about the heating 
mechanism. 

Having specified a radial temperature law, the spectral energy distribution of the star-disk 
system may then be calculated by adding the disk spectrum to the stellar spectrum, where the disk 
spectrum results from a sum, from a minimum radius Rmin = Rhole to a maximum radius Rmax, 
of blackbody annuli at temperatures T(r). 
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Since for a given disk inclination the SED completely specifies the model {Tiau and Ruoie)^ 
the corresponding visibility curves may be computed and compared to the data. The result is 
that the models of HSVK make the prediction that these sources should appear unresolved or very 
nearly unresolved to the IOTA. This conclusion may be qualitatively understood as follows: the 
temperatures at the inner edge of the disk that result from the parameters derived by HSVK are 
in the range T{Rhoie) = 1500 — 3100K (2000 K average); and recalling that the peak of emission 
for a blackbody at 2000 K is Xmax = 0.290 x 10'^/2000K = 1.45 firn, it follows that the region from 
which the near-IR emission arises in those models approximately coincides with the inner edge of 
the disk. Using the average inner hole size (0.3 AU) and distance (617 pc) in their sample implies 
an angular radius of only 0.5 mas for the near-IR emitting region. 

However, the inner hole sizes derived by HSVK are under-estimates, due to their assumption 
that the disks are seen face-on. Thus, we must consider whether a solution that matches the 
visibility data exists for accretion disk models with non-zero inclination. For every value of the 
disk inclination, there will in general exist values of Rhole and Tiau such that the SED is well fit, 
since it is always possible to recover the flux lost to the lower apparent surface area by increasing the 
disk temperature, and to recover the characteristic shape of the SED in the near-IR by increasing 
the size of the inner hole. This degeneracy is illustrated for AB Aur in Figure 7, which shows 
equivalent fits to the SED for inclination values of i = 0° and i = 45° (angle between the disk 
normal and the line of sight). Having found Rhole and Tiau such that the SED is well fit for every 
disk inclination between 0° and 90° , our approach is to determine the smallest inclination for which 
there exists a disk position angle such that the visibility data is also well fit. Note that although a 
non-zero inclination results in larger sources of elliptical shape which can reproduce the amount of 
resolution observed by aligning the disk semi-major axis with the IOTA baseline, it also predicts 
visibilities that vary with baseline position angle, a feature not observed in our data set. 

As before, we require that the flux in the disk component matches the observed excess flux. 
The normalized visibility is given by 



The visibility function for the disk component is calculated by summing the visibility functions 
of the successive annuli, which in turn are given by the difference between the UD visibility functions 
corresponding to the outer {vout) and inner (rj„) annulus radii 



V„{u,v) = 



with: Vr'(^,^^) = FT[4"-'=(C,r?)] 



(18) 



'max 




(19) 



R, 



'max 




(20) 
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For a circularly symmetric UD intensity that is thus tilted by angle i and rotated so that its 

semi-major axis has a position angle ^ (measured from North toward East), the visibility function 
is obtained from the familiar face-on result by a coordinate transformation and use of the scaling 
property of the Fourier transform, so that each disk visibility in Equation 20 is given by 



V^^(«,^;;r) = ^[B^mn) + B^Tout)] ■ ^ ■ cos(i) ■ 2 ■ 



with 



f3 = s • • Wcos2(5' - $) + cos2(i)sin^(5' - $) 



(21) 
(22) 



where the angular diameter of each UD is 9 = 2r /d. 

The results of fitting the SED and visibility data to this model are summarized in Table 9 
and Figure 8. In addition to the quantities already defined, the table lists in the last column the 
temperatures at the inner edge of the disk that result from Equation 17. For those cases for which 
a good fit to the visibility data was found, the best-fit parameters are listed in the table and that 
is the solution plotted in the figure. Note however that for AB Aur and MWC 1080-A, although 
there is no accretion disk model which fits the visibility data well, the data covers a relatively large 
range of baseline position angles (45 and 65° respectively) and we have plotted the closest fit that 
can be obtained in order to illustrate the fact that these solutions are rejected on the basis that 
the data does not show the variation with baseline position angle that is expected from the source 
assymetry. For two other sources (MWC 863-A and MWC 614) the model can not reproduce the 
observed visibilities for any i < 90°, but as an illustration, the table parameters and visibility 
curve plotted correspond to a disk inclination of 30°, the expected average for a randomly selected 
sample. Finally, we note that for V594 Gas, as was the case of the ring model, due to the large 
near-IR excess the accretion disk model can not reproduce the observation that the source appears 
unresolved. This source is the exception, in that the closest fit found corresponds to orienting a 
highly tilted disk with its small axis along the projected baseline. 

In summary, of 11 candidates, 6 can be fit with the disk model. Of those, 2 are unresolved 
(T Ori and MWC 147-A) and consistent with a face-on disk. Of the 4 resolved sources however, 
3 require extreme values of the disk inclination (i > 80°). Based on this statistically unlikely 
requirement, we consider the accretion disk a less attractive explanation for the observed emission. 

For the 4 resolved stars for which the visibility data could be fit, we have computed in Table 10 
the temperature at 1 AU that would result from stellar heating only (Tj^^), by setting M = in 
Equation 16, and where we have used our adopted stellar temperatures, the stellar radii derived 
from the solid angle fits to the SEDs (Table 4) and adopted distances (Table 3.2). Comparing 
with the values of Tiau required to match the SED and visibility data, it can be seen that all 
objects except MWC 297 require an additional source of heating. Under the assumption that this 
additional heating mechanism is viscous accretion, the table shows the value of Mace that would 



-36- 



be required, again using Equation 16. The results reproduce those of HSVK for MWC 297, since 
that is the only object which could be fit with a facc-on disk (also, non-black grains are required 
for this object, since wc obtained TijUj < T^^jj). For the other sources, our implied accretion rates 
are higher than predicted by HSVK, by factors of 30.1, 11.5 and 9.1 for V380 Ori-A, MWC 275 
and V1685 Cyg respectively. 




Fig. 7. — Example of fitting the SED data to the point source plus accretion disk model for two 
different values of the disk inclination, i = 0° and i = 45°. 
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Table 9. Accretion disk model parameters I. 



Name 


i 


Rhole 


TiAU 




Good Fit 


FitRMS 


T{Rhole) 




n 


(AU) 


(K) 


n 


to Visibility Data? 


(%) 


(K) 


V594 Cas 


8L) 


1.5U 


2131 


-bU 


JNo 


8.1 


1572 


AB Aur 


82^" 


0.25 


790 


40 


No 


10.9 


2234 


T Ori 





0.17 


540 


any 


Yes'^ 


0.5 


2040 


V380 Ori-A 


80 


0.41 


1300 


30 


Yes 


2.5 


2537 


MWC 147-A 





0.30 


910 


any 


Yes'^ 


4.1 


2245 


MWC 863-A 


30"^ 


0.07 


280 


30 


No 


20.1 


2057 


MWC 275 


85 


0.24 


720 





Yes 


2.1 


2100 


MWC 297 





0.90 


1830 


any 


Yes 


2.5 


1980 


MWC 614 


30^ 


0.12 


380 





No 


40.2 


1864 


V1685 Cyg 


85 


1.05 


2500 


20 


Yes 


2.2 


2410 


MWC 1080-A 


80^ 


2.70 


5350 


50 


No 


8.5 


2539 



^No good fit exists, value of inclination chosen as an illustration as that which 
predicts approximately the correct size, but is ruled out due to resulting asymmetry, 

see text 

'^No good fit exists for any value of the inclination, value of inclination chosen as an 
illustration corresponds to expected average for randomly selected sample, see text 

^Unresolved 



Table 10. Accretion disk model parameters II. 



Name 










Mace 




(K) 


(Re) 


(K) 


(Mq) 




V380 Ori-A 


12190 


1.7 


229 


3.6 


15.1 


MWC 275 


9332 


2.2 


213 


3.3 


1.5 


MWC 297 


33884 


8.8 


2188 


26.5 


a 


V1685 Cyg 


21878 


4.9 


910 


12.6 


57.9 



References. — Stellar masses from Hillenbrand et al. (1992). 
^No accretion required 
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Fig. 8. — Visibility data and point source plus accretion disk models. Open symbols and solid lines: 
long baseline data and models. Filled symbols and dashed lines: short baseline data and models. 
The models are calculated for the parameters of Table 9. The data and models are plotted as a 
function of source hour angle (H.A.), which determines the orientation of the projected baseline. 
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5.2.4- Infrared Companion 

The potential of an infrared companion for reproducing the observed SEDs of HAEBE stars 
was recognized by Hartmann et al. (1993). In this scenario, the companion is embedded in a 
dust envelope and appears as an infrared source. Although the only star for which a clear binary 
detection has been made is MWC 361-A, we must consider the possibility that the other stars appear 
resolved because we have sampled only a few points of what is really the sinusoidal visibility curve 
of a binary system^. 

In general, the normalized visibility of the central fringe (zero OPD) for a binary system where 
the companion is at offset (^2,^/2) is given by 

V^{u, v) = Al + Bl- F2(Au;) + 2-A^-B^- F(Au;) • cos[27r«2 + vr]2)\ (23) 

where the "fringe washing function" F(Au') = sin(7r Az/Aiu) / (tt Az/Aiu) accounts for the reduction in 
amplitude due to, for a given separation in delay Au;(/xm) = Ao(— n^2+i'?72), the finite bandwidth of 
the spectral filter (represented here by a square function of wavcmimbcr width Az/). The amplitude 
parameters in the above equation are the product of the fraction of the total flux in each component 
times their normalized visibility as individual sources 



A. = p ■ K,i (24) 

B, = • ^-.^ (25) 

^v,l + -t'i/,2 

Therefore, there are four parameters in the model: ^2,112, -^u and Bi,, and unless the data 
provides good coverage in baseline coordinates, further assumptions are needed in order to properly 
constrain it. Thus, we have divided our fits to binary models into three categories. 



Binary model I. Projected separations 

In the case that the data are obtained under essentially a constant baseline position angle, at 
each epoch, we find the angular separation 9 between the binary components "projected" in the 
direction of the baseline vector defined by (u^2 + vr)2) = s ■ & ■ cos($ — T) = s ■ 9; where is the 



^Corporon & Lagrange (1999) also surveyed 42 HAEBE stars spectroscopically and detected 14 companions; 7 
through Li I detection and 6 through radial velocity variations. For those stars which are adaptive optics binaries, 
the large separations measured imply that those companions (the "B" components) are unlikely to be responsible for 
the radial velocity variations detected, and that a third star may be involved instead. Since we measure the fringe 
corresponding to the "A" adaptive optics components only, it is plausible that it is the spectroscopic companions 
which are responsible for the low visibilities we observe. 
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true angular separation and T is the position angle of the binary vector (measured from North, 
toward East). Furthermore, we assume that each star is unresolved {V12 = 1-0), and that all of 
the near-IR excess is due to the companion star, so that the amplitude of the sinusoidal visibility 
is given by the flux ratio, F^^2/F^,i = l^t^ the binary components. With these assumptions, 
there is one parameter left to fit, namely the projected angular separation 9. These assumptions 
are also used in the case of the unresolved sources in order to determine a limit to this parameter. 
The results are summarized in Table 11 and Figure 9. 

Binary model II. Limits on companion position 

For eight of the targets, the visibility data remains constant within the errors, even as the baseline 
position angle changes by a significant amount. In these cases, the lack of visibility variation 
provides a useful constraint to the range of possible binary solutions. Without a detection of 
the sinusoidal visibility however, the relative fluxes, and amount of resolution achieved in each 
component are undetermined. Thus, here we also assume that the totality of the near-IR excess is 
due to the companion, and that both are unresolved. The two parameters left to fit are then the 
angular offsets .^2 and 772 of the companion star. Because observations made at different epochs are 
typically separated by several months, during which the companion position may have changed, we 
establish the possible binary solutions using the single-epoch set of data with the largest baseline 
position angle coverage. 

The results of exploring the solutions allowed by the data are summarized in Table 12, and 
include, in addition to already defined quantities, the formal best fit for the angular offsets of 
the companion (columns 3 and 4) and the range in binary position angle corresponding to other 
solutions consistent with the data (column 6). The formal best fits are plotted together with the 
data in Figure 10. 



Table 11. Binary model parameters I. 



Name 


{F2/F^)h-K 


e 

(mas) 


Fit RMS 

(%) 


N 


V594 Cas 


3.3 ± 0.9 


< 0.7 




6 


T Ori 


9.0 ± 3.2 


< 1.7 




2 


MWC 147-A 


1.7 ± 0.5; 5.6 ± 1.2 


< 0.9 




12 


MWC 863-A 


0.9 ± 0.5; 2.7 ± 0.9 


3.7 ± 0.3 


7.1 


18 


V1685 Cyg 


2.4 ± 0.5 


1.96 ± 0.1 


2.8 


3 
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Fig. 9. — Binary model I: Projected binary separations. The dashed hnes axe upper hmits for the 
unresolved sources. The models are calculated for the best fit parameters of Table 11. 





Table 12. 


Binary model parameters II. 














Best Fit 








Name 


{F2/Fi)h;K 




r/2 Fit RMS 


N 


T 






(mas) 


(mas) 


(%) 




(°) 


AB Aur 


2.2 ± 0.8; 4.7 ± 1.9 


1.3 ± 0.1 


6.3 ± 0.04 


3.4 


21 


11.6 


V380 Ori-A 


5.6 ± 1.6 


-8.4 ± 0.4 


15.9 ± 0.4 


1.8 


7 


[-27.8, +12.9] 


MWC 166-A 


0.1 ± 0.2 


-29.9 ± 4.1 


89.7 ± 1.1 


1.4 


8 


[-35.0, +21.5] 


MWC 275 


1.04 ± 0.6 


-0.6 ± 0.2 


3.2 ± 0.1 


1.3 


9 


[-13.2,-11.9] 


MWC 297 


1.3 ± 0.3 


1.6 ± 0.9 


45.8 ± 0.3 


1.9 


10 


[-1.3 ,+8.1 ] 


MWC 614 


0.8 ± 0.4 


4.6 ± 2.6 


45.9 ± 0.4 


2.2 


9 


[-9.9 ,+17.9] 


V1295 Aql 


0.8 ± 0.4 


1.9 ± 1.9 


28.9 ± 0.3 


4.5 


4 


[-7.4 ,+26.2] 


MWC 1080-A 


3.8 ± 0.1; 10.1 ± 3.5 


2.1 ± 0.2 


2.4 ± 0.1 


3.5 


12 


41.2 
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Fig. 10. — Binary model II: Companion offsets. The models are calculated for the best fit parame- 
ters of Table 12. The data and models are plotted as a function of source hour angle (H.A.), which 
determines the orientation of the projected baseline. In these fits only visibility data corresponding 
a single epoch with the highest range of baseline rotation, which provides maximum sensitivity to 
binary detection, is used. 
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In summary of the last two sub-sections, a binary model can be ruled out for AB Aur, since 
there is no solution that reproduces the constant visibility observed on this source over the 45° range 
of baseline rotation. For the other sources, a binary model can potentially explain the visibility 
data. For the case II sources, the fit residuals are somewhat smaller than in the case of the Gaussian 
and ring models, as one would expect since there are two free parameters in this model. However, 
the difference is within the expected calibration errors and we do not believe it to be significant. 
We note that for MWC 166-A the solution with smallest angular offsets corresponds to a separation 
of ~ 0.1", thus it is unlikely that such companion would have been un-detected in the adaptive 
optics survey by Corporon (1998). 

Binary model III. Solution for MWC 361- A 

For this source, the visibility data shows the clear signature of a binary system, detected interfer- 
ometrically for the first time. Observations were made at two epochs, June 1998 and September 
1998; and both H and K' band data were obtained in each epoch. It may be seen from the data that, 
because the maximum visibility is less than 100%, a solution where at least one of the components 
is resolved is required. Therefore, we fit the data to the four parameters A^, B^, ^2 and r/2. 

The four sets of data do not have the same power in constraining the model, since they do 
not sample the visibility curve equally well. For this reason, we proceed by finding the best fit 
solution to all four parameters using the best set, June 1998 H-band data. The position for the 
companion thus found is used as a constraint in the fit to the June 1998 K'-band data. Likewise, 
the values of A,y and Bi, may be assumed to be constant in time, and therefore only the new value 
for the companion offsets need to be found from the more limited data set of September 1998. 
This procedure was iterated, starting with every solution that provided an adequate fit to the June 
1998 H-band data set, and converged to the best fit solution summarized in Table 13 and plotted 
in Figure 11. Note that the best solution is found for slightly different companion offsets at each 
epoch. However, the difference is at the 2a level, and therefore only of marginal significance. 

By substituting the values for Ai, and B^, in Equations 24 and 25, and given that the total 
fluxes in the system are known (6.0 Jy at H and 8.1 Jy at K, see Table 5), we have a system of 6 
equations and 8 unknowns, and therefore the full solution for the individual fluxes and visibilities 
is indeterminate. However, in order to gain some insight as to the amount of resolution on each 
component that is implied, we solve the above system of equations under the two limiting cases that 
the primary or secondary component is unresolved. The results are summarized in Table 14. The 
table also shows the Gaussian model FWHM (6*"^) that corresponds to the visibility determined 
for each component. It can be seen from the table that the implied sizes for the dust envelopes 
around each component are comparable to those obtained for the other stars in our sample under 
the Gaussian model of Section 5.2.1. 

Corporon Sz Lagrange (1999) found some evidence of radial velocity variations for the MWC 361 
system in various absorption lines of He I and Mg II, but did not have enough data to set a period. If 
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it is the star responsible for the radial velocity variations that we have detected interferometrically, 
this system is a prime candidate for continued observations, using both techniques, from which a 
full orbital solution and individual stellar masses may be obtained. 




Fig. 11. — Binary model solution for MWC 361- A. 
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Table 13. Binary solution for MWC 361-A. 



Epoch & Band 


A. 




6 


m 


N 


Fit RMS 








(mas) 


(mas) 




(%) 


June 1998 - H 


0.772 ± 0.005 


0.156 ±0.008 


0.7 ±0.2 


18.2 ±0.2 


27 


2.3 


June 1998 - K' 


0.820 ±0.12 


0.094 ±0.019 


a 


a 


9 


3.2 


September 1998 - H 


a 


a 


2.4 ± 1.1 


17.7 ±0.3 


13 


6.3 


September 1998 - K' 


b 


b 


c 


c 


5 


2.8 



^Constrained to June 1998 - H solution 

'^Constrained to June 1998 - K' solution 
^Constrained to September 1998 - H solution 



Table 14. MWC 361 limiting solutions. 







assume: ' 


= 1.0, 


then: 






(Jy) 


Ff 

(Jy) 


(mas) 


pK 

(Jy) 


Ff 

(Jy) 


^2 


(mas) 


4.6 


1.4 


0.7 4.0 


6.6 


1.5 


0.5 


6.8 






assume: V2 


= 1.0, then: 






pH 






pK 


Ff 




ef 


5.06 


0.9 


0.9 1.9 


7.3 


0.8 


0.9 


2.7 
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6. DISCUSSION 

The results of our study have demonstrated that, contrary to previous behef, most HAEBE 
stars are easily resolvable in the near-IR by the current class of ground based interferometers. 
Thus, the question of the precise distribution of CS material in these young stars is one that can 
potentially be solved using this technique. 

Ideally, one would like to find a single physical model which can uniquely explain all the sources. 
However, we know that at least 1 of the 15 sources is a binary system (MWC 361-A), and moreover, 
independent evidence (spectroscopic and photometric) points to the binary nature of many other 
HAEBE stars. Therefore, although the results of fitting the visibility data to binary models are not 
significantly better than the fits to simpler 1-parameter models, our data are certainly consistent 
with such an interpretation, and it is likely that this will turn out to be the correct explanation 
for at least some of the sources in our sample. That said, we now consider whether the ensemble 
of observations can be understood in terms of a single model of the near-IR excess. We find that, 
although our data are too limited in spatial frequency coverage to uniquely determine a model, the 
results for the ensemble of sources have features that suggest the interpretation of the emission as 
arising in spherical envelopes around the stars. 

The accretion disk model fails in most cases, and therefore may be ruled out as a candidate 
for the physical explanation of the IR excess. The problem is that, in general, the near-IR emitting 
region must be located too close to the star in order to reproduce the observed SEDs, resulting in 
predictions for the visibilities that are high compared to the observations. This remains true even if 
the disks are allowed to be inclined to the observer, and with the exceptions of MWC 297 and two 
of the unresolved sources, T Ori and MWC 147-A, the model can only reproduce the observations 
in some cases with extreme values of the inclination angle and specific values of the disk position 
angle. Moreover, due to the larger inner holes and higher accretion rates implied, an additional 
difficulty arises with these solutions in that they exacerbate the difficulties raised by Hartmann 
et al. (1993) concerning the physical nature of the opacity holes. We note that, alternatively to 
tilting the disks, the visibility data may be matched by considering a much steeper temperature law 
in the inner edge of the disk (T oc r~^, similar to the ring model), so that the H and K' emission 
arise essentially from the same radius. Thus, if an accretion disk is the explanation for the near-IR 
excess, its inner structure must be quite different from the generally accepted range of power law 
models. 

It is also interesting to compare our results with the mm-wave aperture synthesis observations 
of Mannings &; Sargent (1997) (hereafter MS), who concluded that the mm-wave circumstellar 
emission from HAEBE stars was consistent with accretion disks. Of the sources in our sample, three 
were also part of their study: AB Aur (unresolved in continuum emission, resolved in molecular 
line emission), MWC 863 (unresolved in continuum, undetected in molecular lines), and MWC 275 
(resolved in both continuum and molecular emission). For the two resolved sources, the authors 
found an elongated structure with inclinations 76° and 58° and position angles 79° and 126°, 
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for AB Aur and MWC 275 respectively. The presence of an accretion disk was inferred based 

on the derived disk sizes and masses, which would result in extinctions toward the stars greater 
than observed if the material were spherically distributed, and on line velocity maps, which reveal 
velocity gradients consistent with orbiting material in Keplerian rotation. As we have seen however, 
the analysis of AB Aur using our accretion disk model indicates that the observed symmetry is 
inconsistent with the inclination derived by MS. For MWC 275, the accretion disk model could 
fit the data well only for inclinations > 85°, again inconsistent with the result of MS. Therefore, 
if the mm-wave data do imply circumstellar disk geometry, it appears that at least in these two 
examples the inner regions (< 1 AU) do not follow a simple extrapolation from the disk structure 
at hundreds of AU from the star. 

It is clear from our analysis that both the Gaussian and ring models provide good fits to the data 
(the only two exceptions being V594 Cas and MWC 297, for which the ring model fails to reproduce 
the measured sizes and excess near-IR fluxes). Moreover, it is a striking feature of our data set 
that, except for the obvious binary, none of the sources shows any departure from a brightness 
distribution that appears circularly symmetric on the sky. This statement is most significant for 
sources which have observations spanning a range of baseline position angles. Although the rotation 
of the baseline is essentially zero in some sources, a more typical value is 30°, with a maximum 
of 65° in the greatest case. Thus, to be consistent with the ensemble of observations, models in 
which the emitting dust is distributed symmetrically around the star are favored. In the case of 
AB Aur, for which the baseline rotation is among the largest (50°), our results are consistent with 
the overall symmetry of circumstellar material at radii > l" found in visible light Hubble Space 
Telescope coronographic observations (Grady et al. 1999). For most sources however, it is also clear 
that the indication of circular symmetry needs to be established more firmly with additional data 
covering a larger baseline position angle range. 

Physically the observed circular symmetry would correspond to a spherical distribution of dust, 
as has been proposed in models by Berrilli et al. (1992), Hartmann et al. (1993), Miroshnichenko 
et al. (1997) and Pezzuto et al. (1997). This interpretation is supported by the additional result 
that, in the context of the Gaussian model, the measured sizes and near-IR fluxes require the 
emission to have finite optical depths, consistent with the fact that the central stars need to be 
optically visible when viewed through these envelopes. The ring model in this scenario could then 
interpreted to represent emission from a thin spherical shell, as has been proposed for AB Aur by 
Butner et al. (1994). 

With respect to the physical implications of this model, we find that the properties of the excess 
are not strongly correlated with those of the underlying star. We show in Figure 12 the location 
of the stellar photospheres for the sources in our sample in the HR diagram, based on our adopted 
effective temperatures (Table 4) and luminosities estimated as = Air aR^T^j^^ which leads to 
similar values as those in HSVK and van den Ancker et al. (1998). In the diagram, the sources 
which were found to be resolved are shown using open symbols proportional to the measured linear 
sizes; and solid symbols are used for the size upper limits for unresolved objects. It can be seen 
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that although there is a tendency for the excess sources with largest sizes and highest brightness to 
be located in the upper left part of the diagram, the dependence of those two quantities with stellar 
luminosity and effective temperature is rather weak. We also find that there are pairs of sources, 
(T Ori, V380 Ori-A) and (MWC 147- A, V1685 Cyg), in which the stars and IR excess are essentially 
identical, but the sources of the IR excesses must differ in size by more than a factor of two. This 
might suggest that different physical mechanisms are responsible for the near-IR emission in these 
cases, and that there is no single phenomenon which scales with the properties of the central star 
or the magnitude of the excess in a simple way. Alternatively, if the same underlying mechanism 
is at work in all cases, then it must have the property that the same IR excess is produced by 
systems with very different physical scales. Clearly, such conclusions place profound constraints on 
the nature and physical distribution of the circumstellar material. 
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Fig. 12. — Our sample in the HR diagram. The values plotted correspond to the central photo- 
spheres, and the size of each symbol is proportional to the measured linear size of the source of 
near-IR excess. Solid symbols proportional to the upper limits are used for unresolved sources, and 
open symbols represent resolved sources (in dashed lines if the size is a lower limit). Different sym- 
bols are used for oj Ori (square) because it has no IR excess and for MWC 361 (triangle) because 
it is a binary system. 
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